How do Most Planets Form? - Constraints on Disk Instability from Direct 

Imaging* 



Markus Janson^'^, Mariangela Bonavita^, Hubert Klahr'^, David Lafreniere"^ 



ABSTRACT 

Core accretion and disk instability have traditionally been regarded as the two competing possible paths of 
planet formation. In recent years, evidence have accumulated in favor of core accretion as the dominant 
mode, at least for close-in planets. However, it might be hypothesized that a significant population of 
wide planets formed by disk instabilities could exist at large separations, forming an invisible majority. 
In previous work, we addressed this issue through a direct imaging survey of B2-A0-type stars, and 
concluded that <30% of such stars form and retain planets and brown dwarfs through disk instability, 
leaving core accretion as the likely dominant mechanism. In this paper, we extend this analysis to FGKM- 
type stars by applying a similar analysis to the Gemini Deep Planet Survey (GDPS) sample. The results 
strengthen the conclusion that substellar companions formed and retained around their parent stars by disk 
instabilities are rare. Specifically, we find that the frequency of such companions is <8% for FGKM-type 
stars under our most conservative assumptions, for an outer disk radius of 300 AU, at 99% confidence. 
Furthermore, we find that the frequency is always <10% at 99% confidence independently of outer disk 
radius, for any radius from 5 to 500 AU. We also simulate migration at a wide range of rates, and find that 
the conclusions hold even if the companions move substantially after formation. Hence, core accretion 
remains the likely dominant formation mechanism for the total planet population, for every type of star 
from M-type through B-type. 

Subject headings: planetary systems — brown dwarfs — stars: massive 



1. Introduction 

Direct imaging searches for exoplanets allow to 
constrain the population of w ide and massive ex 



dwarf desert (e.g. Grether & Lineweaveil2006l) . In ad- 
dition, studies of population distributions have demon- 
strated that most of the few wider-orbit exopl anets that 



we used this fact to study the f requency of planets 



have been directly imaged so far (e.g. iMarois et al. 



|2010t|LagrM£eetalJ|2010|) are consistent with exten 



formed by disk instability (e.g. iBossI 120031) . which 



are expected t o hav e exactly these properties (e.g. 



Rafikovl I2005L l2007h . As we noted in that study. 



the population of close-in exoplanets discovered by 
indirect techniques very likely formed by core ac- 
cretion, given lines of evidence such as the correla- 
tion of plane t frequency with stellar metallicity (e.g. 
Fischer et alj l2002i) and the existence of the brown 
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sions of distributions de fined by the closer-in sample 
( Crepp & JohnsonlbOl 1 ). Hence, unless an additional 
wide population of planets formed by disk instabil- 
ity exists in even larger numbers than the presently 
known population, core accretion is the dominant 
mode planet formation. In our previous study, we 
searched for substellar companions around B-stars 
with direct imaging and compared the results to pre- 
dictions of disk instability formation, in order to test 
the presence or absence of such a second population. 
We found that <30% of the stars form and retain plan- 
ets, brown dwarfs and very low-mass stars by this 
mechanism, at 99% confidence. Meanwhile, the fre- 
quency of planets likely formed by core accretion is 



1 



^34% ( Borucki et al. 2011 ^ which indeed implies 
that this is the dominant formation mode for the total 
planet population. 

Still, some uncer tainties remain. For instance, the 
Janson et al. ( 201 iM study was based on B-stars. Part 



of the reason for this was that the massive disks of 
B-stars should provide a very hospitable environment 
for disk instability formation. On the other hand, if 
for some unknown reason it should be the case that 
B-stars happen to be particularly poor hosts for this 
mechanism, then the comparison to the core accre- 
tion population around lower-mass stars could be ir- 
relevant. Hence, in this new study, we extend our 
previous analysis to include FGKM-stars, by using 
archival data from the Gemini Deep Planet Survey 
(GDPS, iLafreniere etaD l2007h . This data set has 



already been analyzed for providing statistica l con- 



straints on planet distributi ons (e.g. Lafrenie re et al 



20071; iNielsen & Closell20lb() . but not in the context of 



disk instability formation. Furthermore, an additional 
uncertainty in the original study was the outer disk ra- 
dius, which was set to 300 AU based on the typical disk 
sizes in the few cases where they have been directly 
measured (e.g. Pietu et al- 2007.) . Here, we study what 
happens if the actual typical disk size is radically dif- 
ferent, testing sizes all the way from 5 to 500 AU. 

2. Analysis 

Our procedure for determining formation limits 
from disk instabilities and evaluation of detectability 
with respect to the detect i on lim its has already been 
described in IJanson et al.l (1201 Ibl) . here we just sum- 
marize the main concepts, and introduce the features 
that are unique to this study. 

The formation limits are set by the fact that two cri- 
teria need to be simultaneously fulfilled for a planet 
to form by disk instability: The Toomre parameter 
([Toomre 1981) needs to be small enough for an in- 
stability to occur, which translates to a minimum mass 
of the resulting planet (referred to as the 'Toomre cri- 
terion' here), and the cooling timescale needs to be 
short enough for fragmentation to occur (e.g. lGamrni3 
20011), which translates to a minimum semi-major axis 
of the resulting planet (called the 'cooling criterion' 
here). These criteria taken together delimit an allowed 
parameter range in mass versus semi-major axis space, 
within which planets and brown dwarfs can form. In 



1 Probably even >50% iMavor et alj201ll) . 



addition, only companions with masses <100 Mjup 
are considered, and the outer disk radius is set to 
300 AU by default (but we test radii all the way from 
5 to 500 AU, as described below). An additional cri- 
terion that we consider here is an upper mass limit on 
the initial fragment from the local available mass, as- 
suming a total disk mass of 0.5 Af^. This limit was 
in principle already in place in our previous study, 
but did not have an impact in practice since it was 
>100 Mjup in the semi-major axis range where for- 
mation was otherwise allowed. In this study, where 
the stellar masses are lower, it does have some (though 
very minor) practical impact. The formation limits are 
generated individually based on the mass, metallicity, 
and initial luminosity of the star In the majority of 
cases where the considered planets would orbit a sin- 
gle star or one component of a binary, the values of the 
parent star in question are used, but for the 14 cases 
where the planets would be circumbinary, the sum of 
masses and initial luminosities of the two components 
is used. The critical s emi-major axes for binari es are 
determined following iHolman & WiegertI ( Il999l) . For 
the cases of unknown binary orbits, we set the semi- 
major axis to 1.31 times th e projected separation fol- 
lowing lFischer et a l . 72002"), and the ecce ntricity as the 
mean of the iDuquennov & Mavon (1199 ih distribution, 
which is 0.36. These critical semi-major axes count 
effectively as additional formation limits in the cases 
where they apply, i.e. if a binary has a critical semi- 
major axis of 50 AU within which stable orbits are not 
possible, then in order for a simulated planet to be al- 
lowed to form, it must fulfill a > 50 AU, in addition to 
the general formation criteria described above. Such 
limits were applied in the 19 cases where binarity sets 
constraints within 500 AU, the semi-major axis cut- 
offs are summarized in Table [T] 

The term 'initial luminosity' in this context denotes 
the luminosity that the star would have had at the point 
when the hypothetical planet formed. This becomes 
a major uncertainty for stars of Sun-like mass and 
lower, because in the time frame that disk instabili- 
ties could be expected to occur, the star is still in its 
Hayashi contraction phase, which means that the lu- 
minosity evolves rapidly with age. In order to remain 
conservative and evaluate this uncertainty, we gener- 
ate two sets of formation limits for each star, based 
on two extremes of initial luminosity, which we refer 
to as cases 'LI' and 'L2'. The two cases correspond 
to timescales which span an order of magnitude, ba- 
sically corresponding to the full TTauri phase of star 
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formation at ages 0.2 Myr to 2 Myr, and evolution- 
ary models ( Siess et aP 12000 *) are used to infer the 
luminosity from a given mass and age. Here, 'LI' 
corresponds to the earlier formation, and thus the 
higher luminosity, and 'L2' to the later formation. 
Stellar masses are inf erred from the spectr al types 
using the relation of Kraus & HiUenbrandl feoOTh 



The metallicities are the same as in 



Lafreniere et al 



(l2007h in the cases where they are listed there, and 
in some cases ne w metallicities from r ecent literature 
have been added ('Alonso et al.' '1996'; 'Ammons et alj 
2006; Paulson & Yelda 2006; Mishenina et al. 20oJ_ 
Holmberg et al.l2009Hjenldns et alJ2009l:ISoubiran etaP 



2010h . [n the few remaining cases where the metallic- 
ity remains unknown, we have assumed Solar abun- 
dance. 

In order to evaluate detectability within the forma- 
tion limits, we use the MESS code (Bonavita et 



201 In to generate 10000 orbits each in grid points of 
5 AU spacing in semi-major axis and 1 Afjup spacing 
in mass. The eccentricities are uniformly distributed 
between 0.0 and 0.6 (the 'eccentric' case of our previ- 
ous study). The projected separation and mass for each 
given orbit is compared to the detection limit from the 
images, and the fraction of orbits that are detectable at 
each given grid point is calculated. The total detection 
probability for the star is then calculated in two differ- 
ent ways: The 'uniform' case in which the planets are 
assumed to be uniformly distributed across the entire 
allowed formation space, and the conservative 'mini- 
mum' case in which the planets are assumed to be uni- 
formly distributed along the minimum mass curve. In 
both of these cases, the semi-major axis distribution is 
assumed to be uniform. The detection hmits are deter- 



mined in the same way as in iLafreniere et al.l (l2007h 
using the 5 a contrast limit translated into mass with 
CON D-based models dAllard et al.l200ll:lBaraffe et al 



I2OO3') for temperatures below 1700 K and DUSTY- 
based models (Chabrier et al. 2000) for higher temper- 
atures. The largest uncertainty in the detection limit 
determination is the age of the individual stars; indeed, 
for a couple of the stars in the GDPS sample, the age 
uncertainty spans more than an order of magnitude. In 
order to fully account for this uncertainty, we evaluate 
two extreme cases separately: One in which the lower 
age limits are used for all of the stars in the sample 
(referred to as the 'young' case), and one in which the 
upper age limits are used (the 'old' case). Examples of 
formation limits and how they relate to detection limits 
are shown in Fig. [T]and Fig. |2] 



100 

90 
80 

: 70 

60 
50 
40 
30 
20 
10 








/Cooling 


Initial 






/ criterion 


mass 




Detection 


\ >^ 


limit 




limits 


] ^ 








1 Allowed 








' — >■ formation 








range 










-— "^Toomre 








criterion 









50 



100 

Semi-major axis (AU) 



150 



200 



Fig. 1 . — Formation and detection limits for HIP 7235, 
which is a KOV-star with a wide possible age range of 
100-1350 Myr. The formation limits correspond to the 
'LI' case. The two detection limits correspond to the 
'old' (upper) and 'young' (lower) cases, respectively. 
As most stars in the sample, the observations have a 
very high degree of completeness, regardless of age. 
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Fig. 2. — Same as Fig. [T] but for the case of 'L2'. 
The completeness remains high, regardless of initial 
luminosity conditions. 
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The upper limit on the frequency of planets formed 
by disk instabilities, /max, is calcul ated f or the full 



population in the same way as in Janso n et al.l (l201 lb), 
using binomial statistics and Bayes theorem. It is 
found that for an outer disk radius of 300 AU, less 
than /max — 7.7% of the stars form and retain planets, 
brown dwarfs and very low-mass stars, at 99% con- 
fidence, in the 'minimum', 'old', 'LI' case. This is 
the most conservative case, i.e. the one that yields the 
highest /max- The lowest /max at 99% confidence is 
acquired in the 'young', 'L2' case where /max — 6.6% 
for both the 'uniform' and the 'minimum' distribution. 

As a further step, we evaluate /max as a function of 
the outer disk radius rout, for radii of 5 AU to 500 AU 
in steps of 5 AU. In the case where rout is smaller than 
the inner edge of planet formation, the individual de- 
tection probability is simply counted as 1 for the pur- 
pose of calculating /max, since disk instability planets 
can't form under this circumstance in the first place 
(in other words, setting /max to 1 implies that we are 
100% certain that no undetectable gravitational insta- 
bility planets could have formed in these cases). The 
results of this test are shown in Fig. |3]and Fig. |4] It 
can be seen that under any circumstance and for any 
outer disk radius below 500 AU, it is always the case 
that /max < 9.9% at 99% confidence. 

The curves in Fig. |3]and Fig. |4]have four distinct ar- 
eas: For very small disks with outer radii of a few 10s 
of AU, planets simply cannot form by disk instabili- 
ties, so the maximum frequency is low. It then sharply 
peaks around ~40 AU, where planets can form and 
simultaneously remain hidden from detection through 
projection effects in some instances. This is followed 
rapid decrease toward outer disk radii of '--^200 AU 
as a continuously decreasing fraction of the formed 
companions can hide through projection. A slower in- 
crease then follows for yet larger radii out to 500 AU, 
as an increasing fraction of the formed companions fall 
outside of the field of view. 

3. Discussion 

The result of this study s i gnifica ntly strengthens the 
conclusion of iJanson et al.l ( 201 lb ) and extends it to 
FGKM-type stars: Core accretion, rather than disk in- 
stability, is likely the dominant form of planet forma- 
tion. Below we provide a discussion of various aspect 
of this result, including several possible caveats that 
need to be taken into account. 
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Fig. 3. — Maximum frequency /max as function of 
outer disk radius rout (99% confidence), for the 'mini- 
mum' distribution case. All 4 possible combinations of 
'young'/' old' and 'L1/L2' are plotted. Even in the ex- 
treme case of 'LI' and 'old', the frequency never rises 
above 9.9% within an outer disk radius of 500 AU. 
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Fig. 4. — Same as Fig. [51 but for the 'uniform' dis- 
tribution case. Even in the extreme case of 'LI' and 
'old', the frequency never rises above 9.6% within an 
outer disk radius of 500 AU. 
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3.1. Planet frequencies 

The aspect of planet formation that we have mainly 
addressed in this paper is which mechanism consti- 
tutes the dominant mode for the total planet popula- 
tion. This includes every type of planet, i.e. both 
rocky and icy planets as well as gas giants. In a sense, 
this might be regarded as an inherently disfavourable 
comparison towards disk instability, since it a mecha- 
nism that is not expected to form any sub-giant plan- 
ets in the first place, which are now known to be the 
dominant component of the core accretion population. 
However, such issues will always take place when two 
mechnisms are compared which have different pre- 
dicted population distributions. Indeed, disc instabil- 
ity is expected to have favourable conditions to form 
very massive gaseous planets, which is difficult within 
the context of core accretion. It is fully conceivable a 
priori that disc instability could form a larger number 
of gaseous planets at large separations than core ac- 
cretion forms rocky planets at small separations. Al- 
though the fact that this is apparently not the case is 
unsuprising, given the high efficiency of core accre- 
tion revealed by Kepler and radial velocity searches, it 
is the most relevant question that can be put forth with 
regards to the planet population at large, and one that 
has not been quantitatively addressed until now. As a 
side point, it should be noted that claims have recently 
been made that disk instability can form rocky planets, 
at least in the sense that the accumulation of the rocky 
body is kick-started in the core of a disk instability that 
subsequently fails to form a compact gaseous planet 
(lNayakshinll2010h . We do not make any strong judge- 
ment on the relative credibility of this scenario here, 
other than to say that more theoretical work is proba- 
bly required in order to distinguish it from the outcome 
of other mechanisms, in terms of observable quantities 
and distributions. 

In any case, another relevant question that can be 
posed following this discussion is which mechanism 
is dominant for the formation of gas giant planets 
in particular. This question can also largely be ad- 
dressed on the basis of our results. For instance, 
ICumming et al. ( 2008) find a giant planet frequency of 
10.5% in their statistically complete sample, and pre- 
dict a frequency of 17-20% by extrapolation to larger 
separations. Even if we adopt the minimum value of 
10.5%, this is still higher than the upper limit of 9.9% 
of disk instability planets that we find in our mo st con- 
servative case. More recently, Johnson el al. ^ 201 oh 
performed a similar study where they also study gi- 



ant planet frequency as function of stellar mass. They 
find a lower value of 8.5% in their statistically com- 
plete sample for FGK-stars. This should also increase 
if extrapolated to larger separations, but akeady the 
number of 8.5% is higher than our upper limit for the 
vast majority of our explored conditions, and in fact, 
to get an upper limit lower than 8.5% even in our most 
pathological case we only need to relax the confidence 
limit to 98%. Hence, it holds true that core accretion is 
likely the dominant formation mechanism even for the 
restricted case of only gas giant planets, for Sun-like 
stars. The frequency of giant planets also increases 



strong ly as a function of stellar mass in Ijohnson el al. 
(120101) . However, this also means that it decreases to 



lower masses, and for M-stars it is only 3.3%. Hence, 
for the specific case of gas giant formation around M- 
dwarfs, we cannot conclude which mechanism is dom- 
inant to any significant degree of certitude. We note 
that this distinction in stellar mass is not relevant in 
terms of comparing formation scenarios for the case 
of the total planet population. The planet frequency as 
func tion of mass is esse ntially flat in the Kepler sam- 
ple (iBorucki et al.ll201 1) . Although the frequency of 
gas giants decreases with stellar mass, the frequency 
of lower-mass planets actually increases toward lower 
stellar masses in their sample. 

3.2. Migration 

While the calculations up to this point have eval- 
uated planets at the location where they formed, it is 
important to note that planets can migrate from their 
natal position to other places in the disk. However, 
while this mechanism can hide individual objects from 
detection, it cannot hide a dominant population formed 
by disk instability, a s we will demonstrate below. In 



Janson et aP ( 201 Ibl) . we put forth two arguments to 



this effect: Firstly, that objects subject to type II migra- 
tion are expected to grow rapidly in mass and would 
quickly hit the hydrogen burning limit and become 
stars, and secondly, that the migration rate would have 
to be extremely fine-tuned in order to hide such com- 
panions from both imaging and radial velocity surveys. 
Here, we quantify these arguments further 

The planet-gas interaction that leads to type II mi- 
gration in gaseous disks is expected to cause the planet 
to grow in mass. This is the mechanism by which 
iMordasini et al.l (l2009h explain the existence of the 
majority of massive gas planets in a core accretion 
framework. These authors find a growth that is rapid 
and has a constant slope in logarithmic mass versus 
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logarithmic semi-major axis, such that for every fac- 
tor 2 that the planet travels in semi-major axis, it gains 
a factor 8.8 in mass. Hence, e.g. a 10 Afjup planet 
that migrates inwards from 50 AU will have a mass of 
88 Afjup - above the hydrogen burning limit - by the 
time it reaches 25 AU. The increase in mass makes the 
object brighter, which makes them yet easier to detect, 
and furthermore, quickly brings them into the stellar 
regime where they are considered beyond the range of 
interest for our study. 

Even if we ignore the growth issue, or assume that 
a mechanism exists that can cause planets to migrate 
without gaining mass, the migration rate (i.e. the dis- 
tance within a gaseous disk that a planet can travel 
within the disk's life time) has to be highly fine-tuned 
in order to hide the planet from detection. If the plan- 
ets travel at a low migration rate, they will be unable 
to reach the inner boundary of the direct imaging de- 
tection hmits, so such an event can be excuded on the 
basis of our non-detections. If, on the other hand, they 
travel at high migration rate, then they would end up in 
the inner regions of the system and be detectable by ra- 
dial velocity surveys (or plunge into the star). Hence, 
in this event, a dominant population of disk instabil- 
ity planets would constitute or strongly contaminate 
the observed radial velocity population. As we have 
discussed in Sect. [T] this would be inconsistent with 
the multiple lines of evidence suggesting that core ac- 
cretion is the dominant mode of formation for the in- 
ner population, so such an event can also be excluded. 
This will be quantified in the remainder of this section. 

The issue of exactly how planets migrate in the 
disk is hard to predict, as dif ferences between re- 



cent modelling effor ts show (e.g. Baruteau et al.ll201 it 
Michael et al.l l201 lb . For instance, Michael et al 



(I2OI 1 ) simulate the first few thousand years of mi 



gration after a giant planet has formed. In one of their 
simulations, the planet migrates rapidly inwards by 
6 AU in ^ 10"^ years and then halts, at least momen- 
tarily. If we interpret this in the sense that planets 
typically migrate ^^6 AU and then stop, then this mo- 
tion is too small to hide any substantial number of 
objects inside the inner detection limit of imaging. If, 
on the other hand, we interpret it in the sense that 6 AU 
per 1000 years is the typical migration rate of planets, 
then in a typical disk lifetime of order 10^ years the 
distance traveled would be thousands of AU, such that 
all disk instability planets would be either hot Jupiters 
or consumed by their parent star. 

Given the uncertainties in quantifying the migra- 



tion rate, we have performed an experiment to demon- 
strate that irrespective of this rate, it is exceedingly 
difficult to hide a population of disk instability com- 
panions through migration, given the highly confined 
parameter space that is available for this purpose. We 
do this by performing a series of similar simulations 
to our main simulation, but instead of evaluating the 
detectability of a given planet at its natal semi-major 
axis, we evaluate it at a smaller semi-major axis, given 
by a particular migration rate. For each simulation, 
we choose different migration rates, from 5 AU to 495 
AU over the disk's lifetime, in steps of 5 AU (i.e., 99 
simulations in total). In order to save computational 
resource, we only evaluate detectability in the 'mini- 
mum', 'L2', 'old' case, which is the most conservative 
set of assumptions possible (i.e., all other sets would 
yield tighter constraints on the resulting upper limit on 
companion frequency). If any companion falls within 
the innermost 5 AU after migration (including AU), 
it is counted as detectable by radial velocity surveys. 
In the cases of circumbinary planets, we halt the mi- 
gration at the inner stability boundary, inside of which 
there should be no disk material to migrate in. For each 
simulation, the upper limit on disk instability compan- 
ion frequency /max is evaluated for the full range of 
disk outer radii from 5 to 500 AU, as in our main sim- 
ulation. 

As a result, we have /max as a function of both disk 
outer radius and migration rate in 99 x 100 grid points. 
In 99.3% of these points, /max < 10% at 99% confi- 
dence. In only one single case is /max > 34%, which 
is our reference lower limit for core accretion. This 
is for the case of a disk radius of ~30 AU and a mi- 
gration rate of ~20 AU per disk lifetime, in which 
a fairly large fraction of the objects formed can hide 
in the most elusive range between radial velocity and 
imaging detection. If we relax the confidence limit to 
97%, then /max = 31%, i.e. under 34% also in this 
case. Hence, we can conclude that very fine-tuned val- 
ues for both the migration rate and the disk radius are 
simultaneously required in order to hide any signifi- 
cant number of targets, and even in that pathological 
case, we can still rule out a dominant disk instability 
population to a high degree of confidence. We also re- 
iterate that this is for the 'minimum', 'L2', 'old' case, 
and that the upper limits would be even tighter for all 
the other cases. Finally, we can also note that if we 
e.g. relax the confidence limits to 95%, then 99.6% of 
the grid points give /max < 8.5%, which is the lower 
frequency limit for gas giant planets around FGK-type 
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stars found by iJohnson el al. I (1201 oh as we discussed 
in the previous section, so we can also exclude disk 
instability as a dominant mechanism for the particu- 
lar population of gas giant planets to a high degree of 
confidence, even in a migration framework. 

3.3. Other possible caveats 

Of course, given the complex and rather poorly un- 
derstood nature of planet formation, there are other 
possible caveats that need to be considered. One fre- 
quent concern regards the uncertainty in the mass- 
luminosity models of young planets, in particular 
with regards to ho t- versus cold-start models (e.g. 
Fortnev et al We discussed this i ssue already in 

our previous study (IJanson et al.ll201 lb ), but to reiter- 
ate, cold-start models are constructed to apply specifi- 
cally to the initial conditions in the core accretion sce- 
nario, and do not apply to disk instability scenarios. 
Since we are concerned with the detection of objects 
formed by disk instability specifically here, the cold- 
start models simply have no relevance. Apart from 
this point, we also showed in Jjanson et al. Il2011b ) that 
even planets that presumably did form by core accre- 
tion do not fit the predicions by th e cold-start models. 
Recen t modelling efforts by e.g. ISpiegel & Burrows 
( 201 11) may be able to mitigate this latter point. 

Another possibly relevant issue concerns uncertain- 
ties in the disk instability formation conditions. Our 
'uniform' versus 'minimum' cases largely encapsu- 
late this uncertainty, but one might imagine that the 
detailed physics of the formation process could al- 
ter the bounda ries of the formation space. Indeed, 
recent work by Kratter & Murrav-Clavl (201 iS, using 



a different approach for the cooling and fragmenta- 
tion treatment, implies that formation might occur 
down to slightly smaller separations than previously 
believed. Furthermore, Meru & Bate (2011) have re- 
cently shown in simulations that a change in disk sur- 
face mass density profiles may affect the cooling cri- 
terion, which would also lead to formation at smaller 
semi-major axes. However, we note that even if the 
detailed boundaries are subject to uncertainties, it re- 
mains the case that there is an enormous parameter 
space further out around every star where such forma- 
tion should be able to occur even more easily, yet it 
remains empty. Hence, from a population perspective, 
this issue is unlikely to matter significantly. 



3.4. Particular objects and classes 

Although disk instability does not seem to con- 
tribute substantially to the formation and retention 
of substellar companions, it is entirely plausible that 
individual objects may form in such a way, at the 
frequency level of a few percent or less. Proba- 
bly, the best candidates for such companions, if they 
exist, are the ^-^20^0 Mjup objects that have been 
found at ^50-100 AU around Sun-like stars, such as 
GO Lup B (Neu hause r et al. 2005; Janson et al. 20o5 
McElwai n 61^111120071) or GJ 758 B (iThalmann et alJ 
i2009c ,Currie et al.ll2010l: Ijanson et al.ll201 lah . These 
could have formed in-situ by such a mechanism, 
whereas core accretion can most likely be excluded 
given the large masses (however, regular binary 
formati o n also rema ins a viable option, see e.g. 
Bate & Bonnelll 2005 ). Other interesting cases in this 



regard are those of very wide (>200 AU) and small 
planet-to-star mass ratio (^1%) systems in young s tel- 
lar as s ociations like Upper S corpius dUafreniere et al.l 



20081: [Ireland et all 12011b iLafreniere et al.1 1201 Ih . 



These have too large separations and small masses to 
fulfill the Toomre criterion and form in-situ by grav- 
itational instability. Hence, if we hypothesize that 
they did form by this mechanism, they would have to 
have migrated outwards. If they anyway had to mi- 
grate outwards, then core accretion remains an option 
to be considered as well, although the masses of the 
companions seem to be on the high end of what this 
mechanism can pr ovide. In either case, f or an object 
like HIP 78530 B dUafreniere et al.ll201 ih in particu- 
lar, it is unlikely that it could have migrated through 
interaction with a gaseous disk, since it is located at 
^700 AU and would therefore require a disk size of 
at least that order Rather, it would have had to have 
been kicked out through dynamical interaction with 
another companion in the system. An interesting de- 
tail in this context is that all these very wide and low 
mass ratio systems have been detected in very young 
(~5 Myr) systems, despite the fact that they should 
be detectable also around older stars. This is consis- 
tent with predictions for dynamical interactions by e.g. 



Veras et al. (2009). In those simulations, a population 



of wide companions at hundreds of AU exist at ages 
of a few Myr, which have been kicked out by other 
companions in the system. These wide companions 
are then continuously ejected from the system through 
further interactions, and no longer exist to any signifi- 
cant extent at ages of ^50 Myr and older. 

A related interesting point is the recent possible in- 
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dication of large quantities of planetary-mass objects 
which may reside either at wide orbits or be free- 
floati ng, from gravitational microlensing (Sum i et aP 
201 ih . If real, we note that it is unlikely that such a 
population could have formed by disk instabilities, be- 
cause the estimated characteristic mass of 1 Mj^p is 
lower than the lowest-mass planets (^2-3 Mjup) that 
can form even around M-type stars. The errors corre- 
sponding to a 95% confidence level in the characteris- 
tic mass of the microlensing population do stretch up 
to '-^8 Mjup, which would be a more reasonable match. 
However, with the much larger masses that would si- 
multaneously be able to form through disk instability, 
it is not clear that the resulting mass function would 
be bottom-heavy enough to create the distinct popu- 
lation of microle nsing event timescales observed by 
Sumi et al.l (1201 lb . whereas core accretion would nat- 
urally produce such a feature in the mass spectrum. 
Hence, on balance we do consider core accretion to 
be the more likely origin also for such a population of 
objects. 

4. Conclusions 

We have studied whether disk instability can form 
and retain a dominant population of substellar com- 
panions around FGKM-type stars, through comparison 
of theoretical model predictions to direct imagin g data 



from the 85 stars in the GDPS survey (Lafreniere et al 



20071) . Thi s work builds on a p revious study of BA 



type stars (IJanson et al.l 1201 la) and extends it to a 
wider range of stellar masses. Our results confirm 
and strengthen our previous conclusion, with the find- 
ing that < 10% of the stars can form and retain such 
companions, at 99% confidence, irrespective of outer 
disk radius for radii from 5 to 500 AU. Since the fre- 
quency of planets formed b y core ac cretion is very 
likely > 34% (e.g. Bor ucki etai]|201lb . it follows that 
core accretion is the most likely origin for the majority 
of the total planet population. This conclusion holds 
true at the 98% or higher confidence level even if we 
restrict ourselves exclusively to the formation of gi- 
ant planets, at least for FGK-type stars, where previ- 
ous studies have found minimum giant planet frequen 



ciesof ~9-10% (ICumming et 



am giant p 
"aiu l2008[ 



Johnson el al 



2010 1). In addition, we find that regardless of migration 
rate from 5 AU to 495 AU during the disk lifetime, it 
is not possible to hide a dominant population of disk 
instability objects simultaneously from direct imaging 
and radial velocity observations. 



The Gemini telescope is operated by the Associa- 
tion of Universities for Research in Astronomy, under 
a cooperative agreement with the NSF on behalf of the 
Gemini partnership. M.J. is funded by the Hubble fel- 
lowship. 

Facilities: Gemini: Gillett (NIRI). 
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Table 1 : Binary limits 


GDPS ID 


Other ID 


Acs (AU)" 


acb (AU)^ 


8 


HIP 11072 


5.7 


30.2 


ir 


HIP 13081 


2.4 


26.4 


14 


HIP 14954 


59.1 


391.2 


23 


HIP 30920 


2.2 


15.3 


28 


HIP 44458 


31.2 


238.7 


38 


HIP 54155 


166.3 


1264.2 


41 


HIP 57494 


16.4 


125.3 


48 


HIP 63742 


0.3 


3.1 


55 


HIP 71631 


4.9 


59.9 


58 


HIP 72567 


28.5 


195.2 


59 


HIP 74045 


5.1 


33.7 


61 


HIP 77408 


185.3 


1370.1 


65 


HD 160934 


2.0 


15.5 


67 


HIP 88848 


0.9 


10.5 


68 


HIP 89005 


18.4 


136.1 


78 


HIP 1 1 1449 


85.9 


628.3 


82 


HIP 115147 


135.9 


975.0 


84 


HIP 116384 


7.0 


54.0 


85 


HIP 117410 


11.2 


85.8 


Circumstellar outer semi-major axis 


limit. 




'Circumbinary 


inner semi-major axis 


limit. 





'^The orbital parameters of the binary are from the preliminary spectroscopic fit of lLatham et alj )2002l) . Furthermore, the system is a triple with an 
additional wide companion to this close binary pair This sets a further constraint on the semi-major axis limit, such that circumbinary companions 
around the close pair can only reside between 26.4 AU and 284.6 AU 
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